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I. INTRODUCTION 



As a background independent approach to quantize general relativity (GR), loop quantum gravity (LQG) has been 
widely investigated in the past 25 years jlJ-Q - Recently, this non-perturbatively loop quantization procedure has been 
successfully generalized to the metric f(lZ) theories]!, @, Brsns-Dicke theory Q and scalar-tensor theories [I]- In 
fact, the scheme of these loop quantum modifiedgravity theories can be extended to more general metric theories of 
gravity with well-defined geometrical dynamics [9j. However, to go round the extreme complexity of a full theory of 
quantum gravity, one approach usually taken is to apply the formal quantization prescriptions to symmetry-reduced 
models. These relatively simple toy models could be employed to test the ideas and constructions of the full theory 
and to draw some physical predictions. The so-called loop quantum cosmology (LQC) is such a symmetry-reduced 
model from LQG. We refer to [IM1 for reviews on LQC. Similarly, to further test the constructions and explore the 
physical contents of loop quantum scalar-tensor theories, it is desirable to study their symmetry-reduced models, such 
as cosmological models. Among all scalar-tensor theories of gravity, the most simple one is the so-called Brans-Dickc 
theory which was introduced Brans and Dicke in 1961 to modify GR in accordance with Mach's principle [14j . 

The cosmological models of classical Brans-Dicke theory were first studied in (T^, [l6[ . Then many aspects of Brans- 
Dicke cosmology have been widely investigated in the past decades [13] • The scalar field non- minimally coupled with 
curvature in Brans-Dicke theory is even expected to account for the dark energy problem [l8l - l23| , which has become a 
topical issue in cosmology [24| . It should be noted that the solar system experiments constrain the coupling constant 
uj of the original 4-dimensional Brans-Dicke theory to be a very large number [25l [26| . For simplicity consideration and 
consistency with the solar system experiments, we will only consider the original Brans-Dicke theory with coupling 
constant w ^ — |. 

This paper is organized as follows. The canonical structure and connection dynamics in the spatially flat FRW model 
of classical Brans-Dicke theory is first given in section UU Then we construct the loop quantum Brans-Dicke cosmology 
in section IIII1 where the dynamical difference equation representing cosmological evolution in the quantum theory is 
derived. In section ITVl by simplifying our quantum Hamiltonian constraint, the path integral method is employed to 
obtain an effective Hamiltonian constraint. In the light of this effective Hamiltonian, the effective dynamical equations 
of loop quantum Brans-Dicke cosmology is derived in section |Vj which implies a quantum bounce near to the classical 
big bang singularity. Conclusions and outlooks are given in the last section. 



II. CANONICAL STRUCTURE OF BRANS-DICKE COSMOLOGY 

The original gravitational action of 4-dimensional Brans-Dicke theory reads [13] 

S (9) = TTT^ / d 4 x^[q>K - ^{8^)8^} (2.1) 



16ttG 
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where (f> is a scalar field, 1Z denotes the scalar curvature of spacetime metric <? Mt ,, and uj is the coupling constant. Now 
we consider the spatially flat, homogeneous and isotropic model. According to the cosmological principle, one can 
write the line element of the spacetime metric of our universe as the following standard form, which is the so-called 
Friedman-Robertson- Walker (FRW) metric 

ds 2 = -dt 2 + a 2 {t) (dr 2 + r 2 (dd 2 + sin 2 6d(j) 2 )) 

where a is the scale factor. At classical level, if one assumes that the matter constituent of the universe be some 
perfect fluid, the evolution equations of Brans-Dicke cosmology would read [15[ 

(2.2) 
(2.3) 

and the equation of motion for the scalar field is 

1 d , ; Q\ 8ttG . „ 
- — "TV <M 3 = — — -P + 3P , 
a 6 dt 3 + 2uj 

where a dot over a letter denotes the derivative with respect to the cosmological time t, p and P are respectively the 
energy density and pressure of the fluid. In the case that the matter part is a massless scalar field, because P = p, 
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the above equation will reduce to 



a 3 dt [ 



la 3 ) 



!6nG 



-p. 



(2.4) 



3 + 2oj 

Recall that loop quantum scalar-tensor theories are based on their connection dynamical formalism Q, where the 
phase space consists of canonical pairs of geometrical conjugate variables, SU (2) connection A l a and dcnsitizcd triad 
E b , and scalar conjugate variables (4>, w). The Poisson brackets between the canonical variables read 



{Ai(x),E b k (y)} 



S(x,y), 



where k = 8irG. To mimic the full theory, we can do the following symmetric reduction of the connection formalism 
as in standard LQC. We first introduce an "elemental cell" V on the homogeneous spatial manifold R 3 and restrict 
all integrals to this elemental cell. Then we choose a fiducial Euclidean metric °q a i, on R 3 as well as the orthonormal 
triad and co-triad ("e"; °wj,), such that °q a b — ° w a 0a; t>- For simplicity, we let the elemental cell V be cubic as measured 
by °q a b and denote its volume by V Q . For spatially flat FRW model we have A l a = jK^, where 7 is a nonzero real 
number and K l a is defined in [8j]. Via fixing the degrees of freedom of local gauge and diffeomorphism transformations, 
we finally yield the reduced connection and densitized triad as [lfj 



Al 



E b =pV 'v^'e . 



where c,p are only functions of t. Hence the phase space of cosmological model consists of conjugate pairs (c,p) and 
(</>, 7r). The Poisson brackets between them read 



{C,P} = g7) 

{4>,tt} = 1. 



(2.5) 



Note that the new variables are related to the old ones by \p\ = a 2 V 3 and c = (<j>a + I^^VJ/ . 

The Gaussian and diffeomorphism constraints in the full theory have been solved by the symmetric reduction. 
Hence, in the cosmological model we only need to treat the remaining Hamiltonian constraint. Its expression in the 
full theory reads 0] 




+ 2^^ 
= 0. 



-y/qDaD^ 
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In the cosmological model which we are considering, the above Hamiltonian constraint reduces to 



H = - 



"f 2 K(f 



3qp 



= 0. 



(2-' 



(2.7) 



(3 + 2w)0|pp 

Recall that in the cosmological model of GR minimally coupled with a massless scalar field, the scalar field can be 
viewed as an emergent internal time variable. An interesting question arising in our Brans-Dicke cosmology is that 
whether the scalar field nonminimally coupled with the geometry can still be viewed as emergent time. To answer 
this question, we check the evolution equation of the scalar field, 



= {<!>, H} = 



2 k 



(3 + 2w)|p|' 



/ 3cp 

K"f 



•7T0). 



If we could show that p^ — + n<j) is a constant of motion, the scalar field would be a monotonic function with 
respect to the cosmological time. This is indeed the case since 



p<p = H} 



6w5Vbl 



6c7rsgn(p) 



+ 



(3 + 2w)k 7 ¥ (3 + 2w) 7v /H (3 + 2w)|p| 



1 



-H w 0, 
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where sgn(p) is the sign function for p. Therefore we conclude that, although the scalar field is nonminimally coupled 
with geometry in Brans-Dicke cosmology, it can still be viewed as an emergent time variable. 



III. LOOP QUANTIZATION OF BRANS-DICKE COSMOLOGY 

To quantize the cosmological model, we first need to construct the quantum kinematics of Brans-Dicke cosmology 
by mimicking the loop quantum scalar-tensor theory. This is the so-called polymer-like quantization. The kinematical 
Hilbert space for the geometry part can be defined as "H|J n := L 2 (R.Bohr, d^n), where Rsohr and dfijj are respectively 
the Bohr compactification of the real line and Haar measure on it [lOj . On the other hand, for convenience we choose 
Schrodinger representation for the scalar field [12J. Thus the kinematical Hilbert space for the scalar field part is 
defined as in usual quantum mechanics, Hy m := L 2 (R,d[i). Hence the whole Hilbert space of Brans-Dicke cosmology 
is a direct product, = %f[ n ® %kin- Now let be the eigenstates of p in the kinematical Hilbert space "H§ n 

such that 



8ttGjH 
PW = — 7, — MlW- 



Then those eigenstates satisfy orthonormal condition 



(MilMj) = <Wj > 



(3.1) 



where is the Kronecker delta function rather than the Dirac distribution. For the convenience of studying 

quantum dynamics, we define new variables 



2V3sgn(p)fi 



-3 



b := /xc, 



where /i = with A = 4\/3iTj£p being a minimum nonzero eigenvalue of the area operator [27]. They also form a 
pair of conjugate variables as 



It turns out that the eigenstates of v also contribute an orthonormal basis in Hf[ n . We denote \<f>, v) as the generalized 
orthonormal basis for the whole Hilbert space H^® . In (b, v) representation, the Hamiltonian constraint (|2.7[) can be 
written as 



H = — 



/3A6 2 



2k 7 2 ^) (3 + 2w)^(A)*|«| 



ihbv 



+ = 0. 



Now we come to the quantum dynamics. As in usual LQC, we start with the Hamiltonian constraint of full theory. 
However, since what we consider here is the homogeneous universe, the last two terms containing spatial derivative 
in the Hamiltonian (|2.6I) can be neglected. Hence we write Eq. (|2.6[) as H = Y^i=i Hi- The quantization of the first 
two terms in H is similar to that in usual LQC [13] . Thus the sum of the first two terms act on a quantum state 
9(v,<t>) ettgf as 

(Hx+H 3 )V(v,<t>) = \ (sin b)A{v) (sin b)*(u,4) 



where 



Hence the final result is 



A(v)^f(i/, <f>) 



V3A 

4«7 2 



\v\\\v + 1| - \v - 



(Hx + H 2 Mv, <!>) = - + 4, 4>) + f (y)*(y, <j>) + /_ {v)*(y - 4, 0)) , 



where 



/+(«) 



/3A 



|« + 3|-|« + 1| |« + 2|, 



I6K7 2 

/_(«) = f+(v - 4), f (v) = -f+(v) - /_(«). 
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Now we turn to H3, H4, H$ terms. Note that here we need to quantize the term K l a Ef . Due to the spatial flatness, 
we have K l a Ef = ^A l a Ef . In the cosmological model, this term can be reduced by 

1 „ 3 Snhbv 
-AlE?~*-cp=—r—. 

7 74 

Because we use polymer representation for geometry, there is no quantum operator corresponding to connection c as 
in standard LQC [12| ■ Hence we have to replace the connection by holonomy to get a well-defined operator. This can 
be achieved by using the classical identity 

lim _5 Jh = lim ^H__2 r . = ~ CTi 

m->o 4/z m->o fj, 

where Tj = — \oi, Oi is Pauli matrices and the holonomy is defined by 

h f) :=C0S ^ 1 + 2s m^T l . (3.2) 



sin(/^c) 

the symmetry- reduced expression of the sum of H3, H4, H$ terms becomes 



here 1 is the 2x2 unit matrix. Thus, according to Eq. (|3.2l) we can replace connection c by holonomy __ki£i Then 



2v3k .Zhbv 2y3t / ,3ft. , sin(6) I v| sin(fe) „ , w 37i N . ... irdn 
(— + irej)) 2 - ■ J —3 (-) 2 ^^J Li + 2 S gn(p)( — ) sin(6)^ + 



(3 + 2w)0(A)5|t;| 4 T/ (3 + 2w)(A)i V 4 ■ •■ |,< 

Let /3 = 3 + 2cj. Based on the above discussion, the action of H3 on a quantum state read 

#3*(<M = 7^3- sm(&)Msin(&)*(«M 

v '" ' '" ' [|u + 2|#(&« + 4)-2|v|*(0,v) + |t>-2|#(>,v-4)] 



2^(A) 

Similarly, for f/4 term we have 



/3(A)* V 4 
Also the action of fig takes the form 



#4*(&») = — — j (^J 2s 5 ra(p) Bin(6)7r*(^») 

2V3k /3ft\ ^ . Nr 9^(0,w + 2) a*(0,u-2). 
~ nr T hs 9 n KP)[ in. 7H. 1 



/o" — — - 
/3(A) 2 



where 



/3(A) ! 

s(«) = (|) 3 H||« + i| 1/3 -|«-i| 1/3 | 3 . (3.3) 

Thus, the Hamiltonian constraint (|2.6[) has been successfully quantized in the cosmological model. The Hamiltonian 
constraint equation of loop quantum Brans-Dicke cosmology reads 



(^ J ff i )*(0,t;) = O. (3.4) 



To study the effective theory of loop quantum Brans-Dicke cosmology, we also want to know the effect of matter 
fields on the dynamical evolution. Hence we now include an extra massless scalar matter field cp into Brans-Dicke 
cosmology. Then classically the total Hamiltonian constraint of the model reads 



H = -^S 2V3. 3^ 2 + ^L =Q 

2k7 2 </> (3 + 2w)0(A)i|u| V 4 Y> \v\{A)r 



6 



where p v is the momentum conjugate to (p. In the quantum theory, the whole Hilbert space now is a direct product 
of three parts, "H^.°* ai = T~L^ in <& ® 'H™ l ^ tter . Here the kinematical Hilbert space for the matter part is also defined 
as in usual quantum mechanics as H™^ tter := L 2 (R, dfi). The action of the Hamiltonian of matter field on a quantum 
state y?) € ?4°* a/ reads 

f m wfu, d>,tp) = T n B(v) — -= . (3.6) 

(A)l ; (A)! <V v ; 



IV. EFFECTIVE HAMILTONIAN OF BRANS-DICKE COSMOLOGY 

To test the robustness of key features of loop quantum cosmology, a simplified soluble model of LQC was proposed 
in Ref.fHI]. In the simplified model, one first adapted the classical theory to the scalar matter field time by the 
following "harmonic gauge" form of spacetime metric: 

ds 2 = -a 6 (r)dT 2 + a 2 (r) (dr 2 + r 2 {d9 2 + sin 2 6W0 2 )) . 

Then the Hamiltonian constraint p.5[) becomes 

A 2 6 2 w 2 2k ,ihbv ,. 2 o „ ,„,s 

In the corresponding quantum theory, we denote quantum state ^(v) = ^S(v,(j),ip) for short. Then the simplified 
Hamiltonian constraint equation reads 

d 2 ^(v) 

{ > -6*(u), (4.2) 



where 



d(p 2 



A 2 

®^( v ) = o w7 v i( v + 2 )*( u + 4 ) ~ 2v9(v) + ( v ~ 2 )^( v ~ 4 ) 



— u [(u + 2)*(v + 4) - 2v^(v) + (w - 2)*(i> - 4)] 



2(3h 2 (f> V 4 

^ ) nv [V(v + 2) - V(v - 2)] + ^rTt(jm^(v) 



iph 2 V 4 y v /J ph 2 

= (YA)*(v). (4.3) 



Thus we get a Klein-Gordon type equation for the quantum dynamics of Brans-Dicke cosmology coupled with a 
massless scalar field. Note that the constraint equation (|4.2I) could also be reduced from the Hamiltonian Q3.4[) and 
p.6[) by the replacements [28j : 

\v\ 

and 

The first replacement amounts to assuming O(-A-) <C 1, which then implies the validity of the second replacement. 

The effective description of LQC is a delicate and topical issue since it may relate the quantum gravity effects 
to low-energy physics. The effective Hamiltonian of LQC are being studied from both canonical perspective [291432T ] 
and path integral perspective [33j- 37] . With the help of the Hamiltonian constraint equation (|4.2p . we now derive an 
effective Hamiltonian within the timeless path integral formalism. In timeless path integral formalism of our model, 
the transition amplitude equals to the physical inner product [33|, |34|, i.e., 

A t i a (vf,$f,(pfi Vi,<pi,tpi) = {vf,$f,<pf\vi,</>i,<pi) p hy = lim / da{vf,(f)f,(pf\e ia \vi,(j)i,(pi), (4.4) 
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where C = 8 + pt,/H 2 . As shown in Refs.[34l. l35j. by multiple group averaging and complete basis inserting, we will 
need to calculate 



{vf,<f>f, <ff\e » =1 \vi,cj)i,(pi) = 



v N _ 1 ,...v 1 



0N-1- 



/ dip N -i...dtpi Y\{ip n \{4>n\{VnW ta,lC \Vn-l)\<l) n—l) \(pn— 1/ • 

(4.5) 



Since the action of the constraint operator C has been separated into gravitational part and matter part, we could 
calculate the exponential on each kinematical space separately. Then, for the matter part one gets 



((fin\e ea "^\(pn-l) = / dp V „(<y5nb y „)(p V „|e" Q "fi^|^„-l) 



2irh 



dp Vn e 



(4.6) 



For the gravity part, we first use the following identity 

M<M(v„|e- l£Q "%„_i)|0 n -i) = -*ea»X) / d<t>n{<t>n\{v n \®i\v n -i)\<i> n -i) + 0(e 2 ). (4.7) 



Then, the matrix elements of <di can be calculated separately by using Eq. (14.3p . We have 



d(/>n{<l>n\(Vn\@l\Vn-l)\<l>r, 



A 2 



8K7 2 S 2 (/)„_i 
1 

2^h 



2 

0« — 1 ■ 



where 7r„ is the momentum conjugate to </>„. Similarly, we can get 



d(j> n (<f>n \{v n I ©2 \v n - 1) |0n-l) 



3ft V 1 



2/?n 2 V 4 7 27rfr 



*n-*n-l >. 1 



— I ^-r / d4> n d-K n e l,i{ - & > '- v n -i 



; «i 1 ,«ti-i+4 



and 



d<t>n (4>n \{Vn I ©3 1 «n-l) |<£n-l 

2k /3ft\ 1 



V 4 / 2ttS 



At last we have 



/3ft 
_K 1_ 



, , rf<^)„ ( / d7r„(0 n ||7r„)(7r„|7r 2 ^|0„_i) + / d7r„(^„|^7r 2 |7r„)(7r n ||0„_i) 



d7r„e 4e "fi l( " «" + 0„_i)7r 2 . 



Taking account of above results and the formula 



i r 



<W„-i+4 ~ 25 Wb ,„ b _ 1 + * OB ,„ n _ 1 _ 4 = - / d6„e- 4b "( u "-""- 1 )/ 2 4sin 2 (&„) 
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Eq. (14.7[) can be expressed as 



M0„IK|e- l£Q ^k„i)|</> n _i) 



1 



1 — ia n e 



A 2 



8n"f 2 h 2 (j), 



V n + V n -i . j 

-v n -i 7, 4 sin b„ 



+ ia n e- 



3ft 

T 1 ''" '" 



V n +V n -1 A . 2/t \ 2ft /3ft \ _ . , v, ft , 



4sin (6„) + ia n£jj^2 7r « l; »-i 2 * sin ( fe ™) + ia " e ^2 (0« + 0n-l)*Vi 



~2/3ft 2 < ? (>„ V 4 y 2 
Collecting all the above ingredients the transition amplitude can be written as 

Ati s (vf,<t>f,(pf, Vi,(t>i,ipi) 



N 



= lim liui | e JJ^ ) / da??... / cfai x / dipN-i---dipi ( ^jj- 



N poo 

■2*1,' I 



dp VN ...dp Vl 



/OC 
dcj>N-i 
-CO 



\27rft 



dbj^...dbi 



7JJV-1,...,«1 



JV 

■ JJ^ exp ze 

n=l 

K 



fV„ <Pn - ¥>n-l . 0n ~ <t>n-l K V n - U n _i 



3ft 



V n -1- 



ft 



a, 



h 2 &K,~f 2 h 2 (f)i 



-V n -1- 



V n + V n -1 



4 sin b„ 



4sin 2 (&„) + - 



Finally, by taking the 'continuum limit' we can get a path integral formulation as 

Ati s (vf,(t>f,(pf, Vi,<f>i,(pi) 

r r r : /•! 

dr 



c Da Dip Vp v / D(j) / 2?7r Vv Vb exp 



■ . , ; hb ( p% 



2 k 



A 2 2 . - 2k /3ft V 2 • 2/^ 4k ^ 3fi \ • <"* j. 2 



where c is an overall constant. Hence, the effective Hamiltonian constraint in the simplified model can be simply read 
as 

A " 2.2, 2k /3ft . ... \ 2 P% 

C ^ f = ~wm v smb+ m[T sm{b)v+n v + f- 

It is easy to see from above expression that the classical Hamiltonian constraint (14.11) can be recovered from C e f / in 
the large scalar limit as sin b — > b. Therefore the above quantum model has correct classical limit. On the other hand, 
if one wants to achieve the effective Hamiltonian constraint for the original model of previous sections, the proper 
time of isotropic observers should be respected. Then the factor A has to be multiplied to C e //. We thus obtain 



H F = — 



/3A 

2"/ 2 K(j) 1 



v \ sm 2 b 



2v/3k /3ft . \ 2 A 3 / 2 \v\ 

— sm(b)v + TT(f>) + n k p, 



where the matter density is defined by 



v% _ §v% 

2| p | 3 ~ w 2 A 3 ' 



2^3 



(4.8) 



Note that the above effective Hamiltonian can also be obtained form the classical Hamiltonian (|3.5[) by the heuristic 
replacement b — > sin b. Hence the classical Hamiltonian constraint can be recovered from the effective Hp in the large 
scalar limit. 



9 



V. EFFECTIVE EQUATION AND QUANTUM BOUNCE 



By employing the effective Hamiltonian H F and symplectic structure of Brans-Dicke cosmology, we can easily get 
equation of motions for v and <fi respectively as 



v = {v,H F } = i-^ | v |sia(b) cos(6) - ^^^ sgn(p) sm(b)v + tt0 j cosiM- 



= {(/>, H F } = 



Ay/in /3fi . 



P&?/ 2 \v\ 



— sm{b)v + ttc, 



(3A 3 / 2 \v\ 



(5.1) 
(5.2) 



Now, let us calculate the evolution oi p^ = ^ sin(b)v + ncj). It reads 

= {P<t>,H F } 

3 A 3 / 2 M 
= -- coa (b)H F +H F --^±±p 

A 3 / 2 \v\ 



2%/3 



Hence p^, is a constant of motion when p = 0. The combination of Eqs. (|5.2p and (|5 .3[) gives 

1 d , • o. 2k 



(5.3) 



which is as same as the classical evolution equation (|2.4[) for the scalar field <j>. However, the evolution equation 
corresponding to Eq. (|2.2|l is modified by the quantum correction, since the combination of equations (|5.1j) and (|5.2p 

gives 



■ sin(b) cos(6) + 



2V3k [Zh 



PAV^v \~ S ' m{b)v + ^ H 1 " cos ^ 



-i 2 



' V jA 2y/A 

On the other hand, the effective Hamiltonian constraint H F = can be rewritten as 

3 sin 2 (b) pj> 2 



-1 2 



(5.4) 



P = 0, 



which gives 



sin 2 (b) = 



(5.5) 



where we defined an effective matter density p e = + 0p and p c = = 32^^—5^ • Note that Eq. (|5.5[) guarantees 

the positivity of p e . Now with the help of Eq. (|5.5p . Eq. (|5.4l) can be expressed as 



a + 24> 



1 [k ~ p e 

-p e (l 

■ > p c 



20 



(1 



(5.6) 



Note that we also have 



6 ( 2kp 2 



v 2 A 3 \ 



v 2 A 3 



P- 



By using Eqs. (|5 . 2[) . (|5 .3[) and (|4.8p . we can show that P? is a constant of motion since 



PI = ^ 
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Therefore, for a contracting universe, p e would monotonically increase while v decreases. Thus, when p e = p c , we 

have cos (b) = ^Jl — y- = 0. Then, from Eq. (l5.1[) . we can easily get v = 0, which implies a quantum bounce happened 

at that point. Hence, just as the LQC case of GR, the big bang singularity of Brans-Dicke cosmology can also be 
avoided by its loop quantization. 

We end up this section with following two remarks on the effective equation (|5.6p . (i) In the special case of = 1, 
we have p e = p and = 0. Then Eq. (|5.6p reduces to the well-known effective Friedman equation of LQC as 

a J 3 p c 

(ii) In the classical limit when p e <C p c , the terms in Eq. (|5.6p can be neglected, and hence it reduces to the 
evolution equation (|2.2p of classical Brans-Dicke cosmology as 



1 k /30 np 




2 i P " 120 2 ' 30' 



VI. CONCLUDING REMARKS 

To summarize the results in previous sections, we first studied the spatially flat FRW model of Brans-Dicke theory 
It turns out that, although the scalar field is non-minimally coupled, it can still be treated as an emergent time variable. 
Hence, in Brans-Dicke cosmology an internal time may come from the gravity rather than an extra matter field. This 
model is then successfully quantized by the nonperturbative loop quantization approach with the Brans-Dicke coupling 
parameter w ^ — |. The Hamiltonian constraint is successfully quantized in this model. Due to the polymer-like 
quantization and the non-vanishing minimal area, the classical differential equation which represents cosmological 
evolution is now replaced by quantum difference equation. In addition, we use the timeless path integral formalism 
and the simplified treatment to derive an effective Hamiltonian of loop quantum Brans-Dicke cosmology. The same 
expression could also be obtained if we took the heuristic replacement of c — > sin (^ c ) in the classical Hamiltonian 
constraint. Hence the quantum theory has correct classical limit. Furthermore, we use this effective Hamiltonian to 
get the effective dynamical equations of the theory, which lay a foundation for the phenomenological investigation to 
possible quantum gravity effects in cosmology. Our analysis indicates that the classical big bang singularity is again 
replaced by a quantum bounce in loop quantum Brans-Dicke cosmology. This result strengthens our confidence that 
the existence of quantum bounce is a universal feature of loop quantum cosmological models. 

There is also interesting situation in loop quantum Brans-Dicke cosmology, which does not exist in the LQC of GR. 
Since the scalar field of Brans-Dicke gravity can play the role of emergent time, there exists a meaningful vacuum 
evolution in loop quantum Brans-Dicke cosmology. In this case, the quantum bounce will still exist even without 

extra matter field. To see this, we let p = 0. Then contributes the effective density p e — &f—, and Eq. (l5.3p implies 
that Pcf) is a constant of motion. Hence Eq. (I5.2p gives oc A and thus p e oc ^ . Therefore the effective density p e 
evolves according to the same rule of the massless scalar matter density p. Then it is easy to see from Eq. (|5.ip that 
a quantum bounce will happen at the point p e — p c . 

It should be noted that there are many aspects of the loop quantum Brans-Dicke cosmology which deserve further 
investigating. For examples, it is still desirable to confirm the effective equations of loop quantum Brans-Dicke 
cosmology from canonical perspective. To confirm the universality of the quantum bounce, we need to generalize our 
scheme to other modified gravity theoris, such as f(R) theories and general scalar-tensor theories. Moreover, since 
our effective equations laid a foundation for the phenomenological investigation to possible quantum gravity effects 
in cosmology, we also would like to further study the cosmological perturbation theory and inflation scenario under 
our framework of loop quantum Brans-Dicke cosmology. We leave all these interesting topics for future study. 
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